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We study the scattering of infrared light by small dielectric core-shell particles taking an AI 2 O 3 
sphere with a CaO core as an example. The extinction efficiency of such a particle shows two intense 
series of resonances attached, respectively, to in-phase and out-of-phase multipolar polarization- 
induced surface charges build-up, respectively, at the core-shell and the shell-vacuum interface. 
Both series, the character of the former may be labelled bonding and the character of the latter 
antibonding, give rise to anomalous scattering. For a given particle radius and filling factor the 
Poynting vector field shows therefore around two wave numbers the complex topology of this type 
of light scattering. Inside the particle the topology depends on the character of the resonance. The 
dissipation of energy inside the particle also reflects the core-shell structure. It depends on the 
resonance and shows strong spatial variations. 

PACS numbers: 42.25.Bs,42.25.Fx,52.27.Lw 


I. INTRODUCTION 

Nanoparticles with a core-shell structure are of inter¬ 
est for various f ield s of applied research. They are rou¬ 
tinely fabrica ted! ^ and have applications ranging from 
nanophotonic^^ to medical treatment.!^ In the field of 
plasmonics,l^for instance, recent attention has been paid 
to nanoparticles havi ng a metallic core and a dielectric 
shelJ^t^ or vice versa.^^^ The reason is the added tun- 
ability of the optical response due to the geometry and 
composition of the particle. One of the most remarkable 
properties of coated nanoparticles is that they may in¬ 
duce transparency in a certain range of wave numbers 
due to the interplay of the dipole radiations of the core 
and the shell.^ii^ 

We are interested in the scattering of light by parti¬ 
cles with a dielectric core and a dielectric coat (shell) 
where both materials show a strong transverse optical 
phonon resonance giving rise to anomalous scattering in 
the infrared occurring for wave numbers where the 
comple x diel ectric function e = e' is" has e' < 0 and 
s" ^ 1 [12IISI In our previous worlP^ we analyzed the ex¬ 
tinction efficiency of this type of core-shell particles with 
an eye on using it in a low-temperature plasma as an elec¬ 
tric probe with an optical read-out. The idea, or iginal ly 
put forward for homogeneous dielectric particlesj ^^ l ^^ l is 
to utilize the blue-shift of the anomalous dipole resonance 
due to the surplus electrons collected from the plasma as 
a diagnostics from which the charge of the particle and 
thus the floating potential at the particle’s position in the 
plasma can be determined. Whereas for homogeneous 
particles the charge-induced shift is most probably too 
small to be of practical importance, core-shell particles 
show a much larger blue shift. In particular the position 
of what we called shell resonance is very charge sensitive. 
Its blue-shift should be measurable for particle radii up 
to 0.5 /rm by infrared attenuation S£ectroscopy as it is 
typically used in plasma diagnostics.^ We also proposed 
to use dielectric particles scattering infrared light in the 
anomalous regime as grains in dusty plasmas and to re¬ 


place traditional force-balance techniqueJ^^^ of deter¬ 
mining the charge of the grains by an optical technique. 

The dielectric core-shell particles suggested for an op¬ 
tical charge measurement contain a core with negative 
and a shell with positive electron affinity. This particu¬ 
lar choice of electron affinities results in a potential well 
which localizes the surplus electrons in the shell of the 
particle. Compared to an homogeneous particle, where 
the surplus electrons are spread out over the whole par¬ 
ticle, the volume charge density of surplus electrons is 
thus enhanced. It is this localization effect which already 
yields larger charge-induced shifts for the anomalous res¬ 
onances of the core. More important however is that the 
core-shell structure leads also to an add itiona l resonance 
not present in an homogeneous particle.l2Sl2ll Within the 
hybridization mode l for the optical response of complex 

nanostructure^22123i additional resonance can be un¬ 
derstood as the antibonding split-off of the anomalous 
resonance of the shell due to its mixing with a cavity 
mode supported by the core. The bonding partner is lo¬ 
cated at smaller wave numbers with the splitting between 
the bonding and the antibonding resonance controlled by 
the filling factor f = b/a which is the core radius b mea¬ 
sured in units of the particle radius a. In our investiga¬ 
tions of charged core-shell particles we found the position 
of the antibonding dipole resonance (which we called shell 
resonance) to be particularly charge-sensitive. 

The purpose of the present work is to analyze the bond¬ 
ing and antibonding resonances showing up in dielectric 
core-shell particles in more detail using a neutral AI2O3 
particle with a CaO core as an example. In addition 
to the results for this particular physical particle we also 
show results for a dissipationless model particle where the 
scattering anomalies can be more clearly identified. We 
demonstrate that the bonding as well as the antibond¬ 
ing resonances scatter light in the anomalous regime. 
This can be deduced from the inverse hierarchy of the 
partial extinction cross sections and the topology of the 
Poyn ting v ector field.^ In contrast to the hybridization 

modeP^J^S] 

we do not use the electrostatic approximation. 
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Figure 1: Illustration of the scattering geometry. The incident 
electromagnetic plane wave is characterized by a Poynting 
vector Sin, an electric held Ein, and a magnetic held Hin. It 
propagates along the 2 -axis with an electric held polarized 
along the x-axis. The index of refraction in the core (i = 1) 
and the shell {i = 2) of the particle is rrii = where Ei is 
the complex dielectric function in region i. The total radius 
of the particle is a while the radius of the core is b. 

Instead we solve the full Mie equation^ 2 lHlIl for the core¬ 
shell particle.l^ We c an thu s gain further insight into the 
hybridization scenaricPHH! by investigating specifically in 
the dipole regime the energy flux inside and outside the 
particle by the methods previously used for homogeneous 
particlesFrom the spatial distribution of the dissi¬ 
pation we can moreover identify regions inside the parti¬ 
cle where most of the energy is deposited which may be 
of technological interest.133 Dissipation depends on the 
character of the resonance-bonding or antibonding-and 
is for both cases highly inhomogeneous, varying by two 
orders of magnitude from one position to another. 

The remaining paper is organized as follows. In the 
next section we recall briefly the Mie theory of light scat¬ 
tering by core-shell particles relegating mathematical de¬ 
tails necessary to fix our notation to an appendix. In 
section Eli we present our results, discussing first the hi¬ 
erarchy of the two series of extinction resonances and 
then-in the dipole regime-the topology of the Poynting 
vector field, the magnitude and orientation of the elec¬ 
tric field, and the dissipation of energy inside the particle. 
Concluding remarks with a summary of our main findings 
are given in section m 


II. THEORETICAL BACKGROUND 

In Fig. we show the geometry of the scattering prob¬ 
lem we consider. An electromagnetic wave, linearly po¬ 
larized in cc-direction and propagating in z-direction, hits 
a core-shell particle located at the center of the coordi¬ 
nate system. The particle is embedded in vacuum. It 
has a total radius a and a core radius b = fa where / 
is the filling factor of the particle. In the limit / —>■ 1 
the particle has thus no shell whereas for / —?► 0 it has no 
core. The core {i = 1) and the shell {i = 2) of the particle 
are made out of different dielectric materials character¬ 
ized by magnetic permeabilities = 1 and dielectric 
functions = e'(A“^) -|- giving rise to re¬ 


fractive indices mi{\~^) = where A“^ = v/c 

is the wave number of the incident infrared radiation, 
V is its frequency, and c is the speed of light. In the 
formulas given below and in the appendix we will use 
w = 27r^' and k = 2 t:X~^ instead of v and A“^, the ab¬ 
breviation ki = krrii, and the size parameters x = 27rA“^6 
and y = 27rA“^a. 

The theory of light scattering by an homogeneous 
spherical particle worked out originally by Mi^^has been 
extended to a core-shell particle by Aden and Kerkei^ a 
long time ago. S ince then the approach has been applied 
extensivelj^l^ and it has also found its way into many 
textbooks.l2SH2Zl Mathematically the Maxwell equations 
for the scattering problem shown in Fig. [^are solved by 
expanding the electromagnetic fields in the three regions 
of interest-the outer space, the shell, and the core-in 
terms of spherical vector harmonics and determining the 
expansion coefficients from the boundary conditions at 
the two interfaces. 

Outside the particle there is the incident and the scat¬ 
tered wave. In terms of vector harmonics the fields of the 
incident wave can be written a^^ 

OO 

Ein = Y.En{M!,ll-tNill) , ( 1 ) 

n—1 

1 OO 

^in = - E (2) 

with expansion coefficients 

= (3) 

n{n + 1 ) 

where Eq is the strength of the incident electric field, 
eventually used as the unit for the field strength, = 
1 is the vacuum permeability, and the superscript (1) 
denotes that the radial dependence of the fields is given 
by the Bessel function j„. The expansions for the fields 
of the scattered wave are 

OO 

A, = E , (4) 

n—1 

j OO 

= Trf E (5) 

n—1 

with a superscript (3) indicating that the Bessel function 
hn describes now the radial dependence of the fields. 

Inside the particle the fields in the shell and the core 


have to be distinguished. Inside the shell, 

OO 

^2 = E - ^9uNili ( 6 ) 

n—1 

+ , (7) 

J OO 

^2 = - E + tfnNili ( 8 ) 

+ , ( 9 ) 
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while inside the core 


III. RESULTS 


Si = ^ - tdnN^li) , (10) 

n—1 

J oo 

= - rTrf E (11) 

W/Tl —' 
n—1 

with the superscripts ( 1 ) and ( 2 ) indicating, respec¬ 
tively, to use the Bessel function and for the ra¬ 
dial dependence of the helds. The expansion coefficients 
o-mbn,Cn,dn, fn, gn,Vn, and Wn have to be determined 
from the boundary conditions at the core-shell and the 
shell-vacuum interface giving rise to a system of eight 
equations for eight unknowns.l^ Once the expansion co¬ 
efficients are known quantities of physical interest can be 
computed. In the following we focus on the extinction 
efficiency, the Poynting vector held, and the dissipation 
of energy inside the particle. 

Of central importance for the characterization of the 
scattering process is the extinction efficiency, that is, the 
scattered and absorbed energy per second divided by the 
incident energy flux per unit area. It is given bjE3 

OO 

Ot = E (12) 

n—1 


with 

Qt"'’ = [( 2 ^^ + 1 )(®" + ^")] ( 13 ) 

the partial extinction efficiency for the multipole reso¬ 
nance of order n where n = 1 denotes the dipole reso¬ 
nance, n = 2 the quadrupole resonance, and so on. 

Whereas Qt is of relevance only in the outer space the 
time-averaged Poynting vector held is of interest in the 
whole space. It can be obtained for the outer and the 
inner regions by inserting the expansions for the helds 
given above into the standard expressioiP^ 


S= 

2 47r 


E X H* 


(14) 


Likewise the dissipation of energy inside the particle can 
be obtained fronPSI 


2 47r 


Lb 


(15) 


using for the electric held and the imaginary part of the 
dielectric function the values applicable to the region un¬ 
der consideration. 

Outside the particle only the scattered and incident 
helds, i?s.in and Hs,in, are present. For the calculation of 
the extinction efficiency and the outside Poynting vector 
held the coefficients a„ and are thus sufficient. In 
order to obtain however the dissipation and the Poynting 
vector held inside the particle the helds Ei ^2 and Hi ^2 
and hence the coefficients Cn,dn, fn, gn,Vn, and Wn are 
also required. In the notation of Bohren and HuffmarP^ 
we list therefore all eight coefficients in the appendix. 


We discuss the scattering of infrared light by neutral 
dielectric core-shell particles with hlling factor / = 0.7 
for two particular cases: (i) a physical Ca 0 /Al 203 par¬ 
ticle and (ii) an idealized model particle. For the former 
we take realistic bulk diel ectric functions for the CaO 
cor^m and the AI 2 O 3 shelP^^^ to be found in our pre- 
vious worlPS and plotted at the bottom of Fig. [ 2 ] while 
for the model particle the dielectric functions of the core 
and the shell are given by 

e' = £00 + ArpQ—L-^ e" = 0 (16) 

with = V jc the wave number of the electromagnetic 
wave, Aij^Q the wave number of the transverse optical 
phonon and Eq and £00 the dielectric constants at low 
and high frequency, respectively. The values taken for 
the core and the shell are listed in Table HI 

In contrast to the dielectric function of the physical 
Ca 0 /Al 203 particle the dielectric function of the model 
particle is real. The model particle has thus no damp¬ 
ing. Since by definition anomalous resonances occur for 
wave numbers for which £' < 0 and e” <C 1 , that is, 
in situations where dissipative losses are much smaller 
than radiative losses,!^ the model particle displays the 
properties of anomalous resonances more clearly than the 
Ca 0 /Al 203 particle which of course has dissipation. In 
particular the disentangling of the complex topology of 
the energy flux is at first easier for the model particle. 


A. Hierarchy of anomalous resonances 

We start with the extinction spectrum in the far-held of 
the particle. In Fig.j^the extinction efficiencies Qt of the 
Ca 0 /Al 203 and the model particle are shown as a func¬ 
tion of wave number A“^ and particle radius a together 
with the dielectric functions of the particles. Here and 
in the hgures below the hlling factor / is always 0.7. In 
the dipole regime, for very small radii, three resonances 
A, B, and C can be identihed for both kinds of particles. 
Resonance C, which we will not discuss in detail, is the 
anomalous dipole resonance of the core while resonance 
B is the one we proposed to use as a charge diagnostics .1^1 
The resonances A and B are of particular interest because 
they belong to a series of resonances lying in the surface 
phonon regime of the shell, that is, in the spectral region 
between AijiQ 2 and A^q 3 , where 3 is the wave 


model 

^TO (cm ^) 

£0 

£00 

core 

300 

3 

2 

shell 

600 

20 

2 


Table I: Parameters for the model particle. 
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Figure 2: (Color online) Top: Extinction efficiency Qt as a 
function of the wave number and the particle radius a 
for filling factor / = 0.7. The anomalous resonance of the 
core is indicated by C while the two resonances appearing 
in the surface phonon regime of the shell are denoted by A 
and B. Bottom: Real and imaginary parts of the dielectric 
functions used for the Ca 0 /Al 2 03 particle and the model 
particle, respectively. 


number of the transverse (longitudinal) optical phonon 
of the shell. The goal of our work is to investigate the 
physical properties of these two series. _ 

From the viewpoint of the hybridization modeP^ESIthe 
two series arise from the mixing of a surface mode of a 
hypothetical particle with radius a made out of shell ma¬ 
terial with a surface mode of a hypothetical cavity of 
radius b made out of core material and embedded into 
a junk of shell material leading in each multipole order 
to a bonding and an antibonding mode. The hybridiza¬ 
tion scenario is buried in the Mie equations but it can 
be made explicit by plotting the difference of the normal 
components of the polarization, P = eo{e — 1 )E, at the 
core-shell and the shell-vacuum interface; Eq = 1 is the di¬ 
electric function of the vacuum. As can be seen in Fig.[^ 
for the particular case of the dipole and quadrupole res¬ 
onances of the model particle, the polarization-induced 
surface charges required to satisfy the boundary condi¬ 
tions of the Maxwell equations at the two interfaces are 
for the bonding resonances in-phase and for the antibond¬ 
ing resonances out-of-phase. The two dipole resonances, 
for instance, can thus be understood in terms of an inner 
and an outer dipole which are oriented either parallel or 
antiparallel. In the vicinity of these two resonances the 
electric field outside the particle is thus the effective field 
of two parallel or antiparallel dipoles. Inside the parti¬ 
cle the field is more complex (see below). Similarly the 
bonding and antibonding quadrupole resonances origi¬ 
nate from quadrupole surface polarization charges at the 
two interfaces which are either in- or out-of-phase. 

The parameters of the dielectric functions of the model 
particle have been chosen in such a way that the series 
of antibonding resonances can be also clearly identified. 
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Figure 3: (Color online) Polarization-induced surface charges 
build-up at the core-shell and the shell-vacuum interface of 
the model particle to satisfy the boundary conditions of the 
Maxwell equations for the dipole (upper two panels) and the 
quadrupole resonances (lower two panels). For the bonding 
resonances the multipolar surface charges of the two interfaces 
are in-phase whereas for the antibonding resonances they are 
out-of-phase. The wave numbers are always slightly above the 
resonance: = 1095.669cm“^ (bonding dipole), A~^ = 

1767.293cm~^ (antibonding dipole), A“^ = 1241.95cm“^ 
(bonding quadrupole), and A“^ = 1682.42 cm“^ (antibonding 
quadrupole). The particle radius a = 0.1pm for the dipole 
and a = 0.6 j-im for the quadrupole resonances with filling 
factor / = 0.7 in both cases. 


Both the bonding and antibonding series stop at the third 
multipole order, that is, at the octupole, as it is also the 
case for the Ca0/Al2 03 particle. But it is conceivable 
to have material combinations and perhaps filling fac¬ 
tors for which the series of sharp resonances continues to 
higher order. The question then arises whether the two 
series of resonances scatter light anomalously. In partic¬ 
ular, it would be interesting to know whether the sharp 
resonances signal the resonant excitation of a particular 
surface mode of the shell and whether the hierarchies of 
extinction cross sections to which the two series give rise 
to are inverted as it would be the case for homogeneous 
particles made out of either the shell or the core material. 

For an homogeneous particle Tribelsky and 
Luk’yanchulP^ have shown that if the size parame¬ 
ter is at the resonance wave number of the order one 
or smaller the main contribution to this resonance 
comes from one resonantly excited surface mode. This 
selectivity remains for a core-shell particle, as we will 
now demonstrate. For that purpose we plot in Fig. for 
the physical and the model particle the extinction effi¬ 
ciencies due to the dipole mode and the quadrupole 
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a) Ca0/Al203 b) model particle 





Figure 4: (Color online) Partial extinction efficiency 
in the vicinity of the bonding (A) and antibonding resonance 
(B) as a function of the wave number and the particle 
radius a. The superscripts (1) and (2) indicate the dipole and 
the quadrupole mode. Higher multipole modes appearing at 
larger radii are weaker and cannot be seen on the scale used 
for the plot. The filling factor / = 0.7. 


mode separately. As can be seen each mode has 
its own distinct set of particle radii and wave numbers 
at which it is resonantly excited. Off-resonance the 
partial multipole extinction efficiencies provide only 
an homogeneous background to the overall extinction 
efficiency. Thus, it is possible to relate to the lowest 
resonance in Fig. an excited dipole mode and to the 
second lowest resonance an excited quadrupole mode. 
For a given particle radius a core-shell particle gives thus 
rise to two wave numbers at which a dipole mode or a 
quadrupole model can be excited. This is in contrast to 
an homogeneous particle where only one wave number 
exists. 

Having identified the multipole order of the resonances 
we can now turn to the extinction cross section Ut = 
TTO^Qt and study its hierarchy. For an homogeneous par¬ 
ticle at shows in the anomalous regime an inverse hierar¬ 
chy, that is, for a given particle radius the cross section 
due to a higher multipol e is la rger than the cross section 
due to a lower multipole. b^b^l 

In Fig. we see for the model particle that the series 
of bonding and antibonding resonances arising in core¬ 
shell particles display the same behavior. At the bond¬ 
ing resonance (left panel) for particle radii smaller than 
a = 1 ^m the cross section at of the dipole resonance at 
« 1050 cm“^ is dominant, at a = 1.3 /rm on the other 
hand, the maximum of the cross section shifted to the 
quadrupole resonance at « 1200cm“^. Yet another 
shift occurred for a = 2.7 /rm, where the maximum of 
the cross section is now at the octupole resonance. This 
inverse hierarchy can also be seen for the antibonding res¬ 
onance (right panel). At the particle radius a = 0.8^m 
the cross section maximum has shifted from the dipole 
mode around A“^ « 1740 cm“^ to the quadrupole mode 
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Figure 5: (Color online) Inverse hierarchy of the extinction 
cross section at in units of for the model particle in the 
vicinity of the bonding (A) and the antibonding resonance 
(B). Horizontal lines indicate radii where the local maximum 
of the cross section shifted from a lower to a higher mode. 
The filling factor / = 0.7. 


around A“^ « 1680cm“^. The next shift of the cross 
section maximum is indicated at a = 1.3/rm where the 
octupole resonance is the most pronounced in the cross 
section. Since the inverse hierarchy is a direct conse- 
quenc e of t he series of resonances in the extinction effi- 
ciency,^i2ll3l Ca 0 /Al 2 03 particle shows the same 

inverted hierarchy of resonances in the cross section (not 
shown in Fig. [^. 

Both the series of bonding and the series of antibond¬ 
ing resonances fulfill therefore the criteria of anoma¬ 
lous light sc atter ing as spelled out by Tribelsky and 
Luk’yanchuk.li21i^ Each series consists of strong narrow 
resonances which are determined by a single resonantly 
excited bonding or antibonding mode and which more¬ 
over give rise to extinction cross sections increasing for 
fixed radius with the order of the multipole. The differ¬ 
ent physical origin of the two series, in-phase vs. out-of- 
phase polarization-induced surface charges at the two in¬ 
terfaces, is reflected in the way the two series shift in wave 
number with increasing multipole order. Whereas the 
bonding resonances shift to higher wave numbers the an¬ 
tibonding resonances shift in the other direction. The an¬ 
tibonding quadrupole resonance is thus excited at lower 
wave numbers than the antibonding dipole resonance. 


B. Topology of the Poynting vector field 

In the previous section we showed that the series 
of bonding and antibonding resonances have anomalous 
properties. We now investigate in the vicinity of the two 
lowest resonances, the bonding and antibonding dipole 
resonances, the energy flux outside and inside the core¬ 
shell particle and compare it with the results Tribelsky 
and Luk’yanchuk obtained for an homogeneous parti¬ 
cle.^ In particular, we consider a model particle with 
radius a = 0.1 ^m and a Ca 0 /Al 203 particle with ra¬ 
dius a = 0.4/rm keeping the filling factor in both cases 
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fixed to / = 0.7. As can be deduced from Fig. |^only the 
dipole resonances are then excited. For the model par¬ 
ticle they are, respectively, at « 1095.669 cm“^ and 
A“Ji « 1767.293 cm“^, where the subscript stands for 
bonding/antibonding (b/ab) dipole (n = 1) resonance, 
while for the Ca 0 /Al 203 particle « 709.63cm“^ 
and « 871.62cm“^. 

Tribelsky and Luk’yanchul^i^ investigated the Poynt- 
ing flux in the vicinity of an anomalous resonance of an 
homogeneous particle and identified its singular points. 
Due to the large changes in the modulus and phase of the 
scattering coefficients in this spectral region the position 
and the occurrence of the singular points is very sensi¬ 
tive to little changes in the wave number of the incident 
light. A singular point is characterized by |5| = 0.1^ For 
the scattering geometry shown in Fig. which is also 
the one used by Tribelsky and Luk’yanchuk,^^ the sin¬ 
gular points are located in the xz-plane. Since in this 
plane S^{x,z) = 0 they are given by the intersection 
of the isoclines Sr{x,z) = 0 and Sg{x,z) = 0. In the 
dipole regime only the coefficients a„, dn, g-n, and 
are resonant. Hence, the orientation of the electric and 
magnetic fields in space is in all three regions determined 
by the vector harmonics and respectively. 

The difference is in the radial dependence, signalled by 
the superscript (...) which may stand for (1), (2), or (3) 
depending on the field under consideration. 

We also identify the type of the singular point which 
in vacuum, where V • S' = 0, can be only a center, a sad¬ 
dle, a saddle-focus or a saddle-node.l^ For that purpose 
we (i) carefully analyze the field lines around the singu¬ 
larity as each type of singularity shows a characteristic 
field distribution and (ii) we check the behavior of the 
magnetic field Hy{x,z) = \H{x, as suggested 

by Schouten.l^ At a center or a focus singularity the am¬ 
plitude of the magnetic field vanishes \H{x, z) \ = 0, while 
a saddle or node point is characterized by 'S/4>(x,z) = 0. 
Singular points due to anomalous resonances all lie in 
the near field since the characteristic spatial scales are 
smaller than the wave length of the radiation.E^ For in¬ 
stance, for the model particle, the largest distance of 
singular points due to the antibonding dipole resonance 
is about r = 1.5 x 10“^ cm while the wavelength of 
the incident light is A = 5.7 x 10“^ cm. Similarly, for 
the bonding dipole resonance, r = 2.5 x 10“^ cm and 
A = 9.13 X 10-4 cm. 


Topology outside the particle 

We start with the outside energy flux. For both the 
model and the Ca 0 /Al 203 particle only the dipole reso¬ 
nances are excited. Since outside the particles the topol¬ 
ogy of the energy flux in the vicinity of the bonding and 
antibonding (dipole) resonances turn out to be qualita¬ 
tively the same, we discuss only results for the latter. 

Figure [^summarizes for the model particle the topol¬ 


ogy of the outside Poynting vector field in the vicinity of 
the antibonding dipole resonance (a = O.l/rm, / = 0.7, 
•^abi ~ 1767.293cm-4). Note, the dramatic modifica¬ 
tions in the topology of the Poynting field due to very 
small changes in the wave number. The top left of the 
figure shows the position of the singular points, labelled 
by Arabic letters, with respect to the particle center. The 
Roman numerals on the top indicate regions with a fixed 
topology. The remaining three panels show the Poynting 
flux outside the particle, symbolized by the two brown 
circles in the center of the plots, for wave numbers rep¬ 
resentative for the indicated topological region. The ar¬ 
rows show the direction of the Poynting vector, while the 
colors correspond to the modulus of the Poynting flux 
normalized to the incident flux of light. In addition the 
Sr = 0 isocline and Sg = 0 isocline are indicated by blue 
and green lines. The singular points are given by the 
intersection of two differently colored isoclines. 

For wave numbers A-^ < A-Jj^, for instance, in re¬ 
gion VIII, the Poynting flux in the vicinity of the particle 
shows a dipole radiation oriented opposite to the incident 
light which is propagating in ^-direction. As an illustra¬ 
tion of this behavior we plotted in the upper right panel 
of Fig. the Poynting flux for A”^ = I767.28cm-4. It 
contains two singular points, one in the upper and one in 
the lower half of the a; 2 :-plane. At the saddle point D in 
the lower half plane (replotted on a magnified scale in the 
inset) it can be nicely seen that the energy flux emitted 
by the particle due to the excited dipoles is oriented oppo¬ 
sitely to the incident energy flux. The field lines emitted 
by the particle (positive 5'r-component) and the field lines 
from the incoming wave (negative S'^-component) meet 
at this point and are then bend over to flow towards plus 
or minus x. In the upper half-plane, that is for z > 0, 
the Poynting flux either passes by or is reabsorbed by the 
particle. This can be seen at point C, which is the Bohren 
saddle point,!^ where the field lines coming from ±.x are 
going behind the particle (positive S'^.-component) or to¬ 
wards it (negative S'^-component). The saddle points C 
and D reflect the strength of the energy flow of the reso¬ 
nant excitation relative to the strength of the energy flow 
of the incoming light. Thus, these points can disappear 
if the dipoles are not emitting enough energy, as it can 
happen if the damping inside the particle is too large. 

On the high energy side of the resonance, for wave 
numbers A-^ > A-Jj^, the dipoles radiate energy parallel 
to the incident light. The Poynting flux thus gets more 
complicated. For example in the topological region III, 
at A-4 = 1767.35cm”^, six singular points, four sad¬ 
dle points A, B, C, D, and two foci points E and F (in 
three dimensions these are saddle-focus points) arise, all 
located in the upper half of the ccz-plane (see lower left 
panel of Fig. [^. There are no singular points anymore for 
z < 0. Since the excited dipoles changed simultaneously 
their directions the 5'r-component of the Poynting vec¬ 
tor is always negative for x < 0 , that is, the energy flux 
along the z-axis is oriented towards the particle. One of 
the Sr = 0 isoclines at which a change of sign in the Sr- 
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Figure 6: (Color online) Top left: Position of the outside singular points in the Poynting field (Arabic letters) of the model 
particle with / = 0.7 and a = 0.1 pm near the antibonding resonance at « 1767.293 cm“^ as a function of wave number. 
For the bonding resonance the results are similar. Spectral regions with a hxed topology are indicated by Roman numbers and 
the shaded area denotes the shell. The other three panels show the Poynting field for three Hxed wave numbers. The color of 
the arrows gives the modulus of the Poynting vector normalized to the flux of incident energy. The intersections of the isoclines 
Sg{r, 9) = 0 (green) and Sr{r, 9) — 0 (blue) determine the singular points. The shell and the core of the particle are indicated 
by brown circles. 


component occurs is located above the particle at z > 0 
and so is now the saddle point D. At this point the field 
lines emitted by the particle and thus having a positive 
Sr-component and the field lines coming from the foci 
points E and F and heading towards the particle owing 
to the negative S'^-component meet and take their way 
towards plus and minus x. This is quite similar to the 
situation at point D in region VIII, since the foci points 


can be regarded as sources of energy in the xz-plane.l^ 
At the Bohren saddle point C the held lines come to¬ 
gether from the left and right focus and then they either 
go behind the particle or towards it, as they did in re¬ 
gion VIII. At z = 0 there are two saddle points A and 
B. To understand their functioning lets take a look along 
the Sg = 0 isocline. For z < 0 held lines which cross the 
Sg = 0 isocline enter the particle while held lines crossing 










































the Sg = 0 isocline for 2 ; > 0 leave the particle. Thus, at 
the singular points A and B field lines leaving or entering 
the particle come close together and form a saddle point. 

Some wave numbers away from in region I, only 

the saddle points A and B remain as can be seen in the 
panel on the lower right of Fig. [^where the Poynting flux 
is plotted for = 1767.38 cm“^. If we would plot sep- 
aratrixes inside the circle defined by the Sg = 0 isocline 
we would find Tribelsky ears.^^ For z < 0 the Poynting 
flux has not changed qualitatively compared to region 
III. But for z > 0 the Sr = 0 isocline disappeared. As a 
result, there are no further singular points other than A 
and B. 

A comparison of our results for the Poynting flux near 
the antibonding dipole resonance of a core-shell particle 
with the results found by Tribelsky and Luk’yanchulsP^ 
near the anomalous resonance of a homogeneous parti¬ 
cle shows that the topologies are the same. The same 
pattern arises also around the bonding dipole resonance. 
Hence, the changes in the anomalous topology of the 
Poynting flux appear for a core-shell particle around two 
wave numbers. By just looking at the energy distribution 
outside the particle one could not distinguish between a 
core-shell and an homogeneous particle. Although in the 
first case two dipoles emit energy, which may have par¬ 
allel or antiparallel orientation, depending on whether 
one is close to the bonding or close to the antibonding 
dipole resonance, whereas in the second case the energy 
is emitted only by one dipole. 


Topology inside the particle 

The core-shell structure of the particle reveals itself in 
the energy flux inside the particle. The reason is the rel¬ 
ative orientation of the inner and the outer dipole. It 
is parallel for the bonding and antiparallel for the anti¬ 
bonding resonance. Below we will verify the orientation 
of the dipoles explicitly by looking at the electric field 
lines. At present it suffices to recall the polarity of the 
surface polarization charges shown in Fig. from which 
the orientation of the dipoles can be also deduced. Unlike 
the Poynting flux outside the particle, which arises from 
the superposition of the incoming field with an effective 
overall dipole field of the particle, inside the core-shell 
particle the energy fluxes of two dipoles which are either 
parallel (bonding) or antiparallel (antibonding) mix. 

As explained at the beginning of this section in the 
dipole regime the singular points of the Poynting flux 
are all located in the a;z-plane. The figures summarizing 
the results for the Poynting flux inside the particle are 
thus constructed in analogy to Fig. We will show the 
position of the singular points with respect to the particle 
center and for representative wave numbers the topology 
of the energy flux. The only difference is that the radial 
distance of the singular points from the center is now less 
than the particle radius a. 


1. Antibonding dipole resonance 

The inside topology of the Poynting flux near the an¬ 
tibonding dipole resonance of the model core-shell par¬ 
ticle is presented in Fig. On the scale of this figure 
the localized surface electromagnetic modes characteris¬ 
tic for the type of light scattering we discuss in this pa¬ 
per can be clearly identified. They are associated with 
the optical whirlpool^^ in the Poynting flux which are 
indicated by light blue circles. The energy is rotating 
around whirlpools and thus Sr = 0. However, since Sg 
has a maximuirP^ at these points whirlpools are not or¬ 
dinary singularities. At the two interfaces the whirlpools 
are oriented in the opposite direction. It is instructive to 
trace the direction of the whirlpools to the orientation of 
the inner and the outer dipole. The polarity of the sur¬ 
face charges at the two interfaces as well as the electric 
field lines to be presented in Sec. |HI C| show that at the 
antibonding resonance the two dipoles are antiparallel. 
Hence, in the shell region the electric field lines belong¬ 
ing to the outer and the inner dipole are parallel. In the 
vicinity of the whirlpools this leads then to a rectifying 
Poynting flux inside the shell, with Poynting field lines 
separated by a S'r = 0 isocline. 

We now go through Fig. step by step. The anti¬ 
bonding dipole resonance is at A“J^ « 1767.293cm“^. A 
representative topology of the Poynting vector field for 
A“^ < A“J;^ is the one of region VHI shown in the up¬ 
per right panel of Fig. for the particular wave number 
A“^ = 1767.289 cm“^. Like before the singular points 
are the intersection points between two lines of different 
color corresponding respectively to the Sr = 0 and Sg = 0 
isocline. Very close to horizontal = 0 isoclines (blue) 
in the negative xz half-plane is also a Sg = 0 isocline 
(green). Both come from the outside and have presently 
no physical consequences. In region VHI there are two 
singular points inside the shell, K and J, and, as discussed 
before (see Fig. [^, two singular points outside the par¬ 
ticle, C and D, from which only D is within the range of 
Fig. [3 The Poynting flux shows that the particle is emit¬ 
ting and reabsorbing light to-and-fro the outside region. 
In the vicinity of the core the Poynting flux in the shell 
resembles the flux around an excited dipole. The core 
is emitting energy into the shell and is also reabsorbing 
energy from it. The singular points inside the shell arise 
therefore from an interference between the electromag¬ 
netic fields of an inner and an outer dipole. 

In region VHI the saddle point K arises from field lines 
emitted from and reabsorbed by the outer dipole meeting 
field lines which are emitted from the inner dipole. The 
field lines take then their way to the left and right, that 
is, to ±x, causing thereby a rectified energy flux. At the 
saddle point J, on the other hand, the rectified energy 
flux is divided, some field lines are reentering the core 
and some are leaving the particle. In the shell region, the 
points J and K operate thus analogously to the points C 
and D in the outside region of the particle (cp. region 
VHI of Fig. 1^. One can hence interpret the reabsorbed 
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Figure 7: (Color online) Top left: Position of the singular points (Arabic letters) occurring in the Poynting field inside the 
model particle (/ = 0.7 and a = O.l/rm) in the vicinity of the antibonding dipole resonance at A“|/j « 1767.293cm“^. The 
different topologies are indicated by Roman numbers and the shaded regions denote the shell. The other three panels show the 
Poynting flux for three wave numbers, each one is representative for a certain topology. The arrows show the direction of the 
Poynting flux at each point, while the modulus of the Poynting vector normalized to the incident flux of energy is encoded by 
the color of the arrows. The intersections of the isoclines S 0 {r, 6 ) — 0 (green) and Sr{r,9) = 0 (blue) determine the singular 
points. The shell and the core are again indicated by brown circles and the optical whirlpools are marked by light blue circles. 


energy flow of the outer dipole as light which is incident 
on the core and interacting there with the inner dipole. 

For wave numbers closer to the resonance the outside 
saddle point D enters the shell region, as can be seen in 
the top left panel of Fig.[^ where for « 1767.29 cm“^ 
the green line corresponding to the singular point D is 
coming up form z/a = —1. In the top left panel of Fig. 
on the other hand, the green line disappears for that 
wave number at r/a = 1. The point D is caused by the 
intersection of the Sr = 0 isocline (blue line) with the 
Se — 0 isocline (green line), which lies on the z-axis. 
With the point D the Sr = 0 isocline is also entering 
the particle. There has to be thus an avoided crossing 


with the 5^ = 0 isocline already existing inside the shell 
because isoclines are solutions of the differential equation 
and a crossing would violate the uniqueness of 
the solution.l^ As a result, see the top left panel of Fig. 
there is only one singular point-the point K-in region 
Vllb, where A“^ < and one singular point-the 

point L-in region Vb, where A“^ > A“J^. The Poynting 
flux for the latter is shown in the lower right panel of 
Fig. [^and will be discussed in a moment. 

The avoided crossing occurs also in region VIIc the 
Poynting flux of it is plotted in the bottom left panel of 
Fig. The Sr = 0 isocline entering the shell from the 
outside goes around it taking its way through the upper 
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Figure 8: (Color online) Position of singular points (Arabic letters) in the Poynting field outside (a) and inside (b) a Ca 0 /Al 203 
particle with / = 0.7 and a = 0.4/rm for wave numbers near the antibonding resonance at j ~ 871.62 cm“^. Roman numbers 
indicate regions with fixed topology which were also found for the model particle. The shaded regions denote the shell. 


half of the xz-plane and causing thereby the saddle point 
K. While the Sr = 0 isocline inside the shell enters the 
core to prevent a crossing and causes thereby the saddle- 
node point M inside the core and the saddle point L inside 
the shell. The saddle point L has the same function as 
the saddle point K. The emitted field lines of the outer 
dipole are reabsorbed in the shell and meet there the 
emitted field lines of the inner dipole. The saddle point 
M, on the other hand, distinguishes between a region 
where both the inner and the outer dipole have already 
changed their directions from a region where they have 
not. In the xz-plane this point is a singularity line. Thus 
to visualize the Poynting flux around this saddle point 
one would have to look at the yz-plane. 

The topology of the Poynting flux in region Vb is 
plotted in the lower right panel of Fig. for = 
1767.2965 cm“^. The singular point L can be clearly 
seen. Since A“^ > both dipoles have changed their 
direction as can be also inferred from the orientation of 
the energy whirlpools at z = 0. For wave numbers well 
above A“J ^ the topology of the energy flux inside the par¬ 
ticle is the one of region IV. There are now two saddle 
points L and J corresponding to the singular points K 
and J in region VIII, respectively. The topology before 
and after the resonance is thus in fact the same only up¬ 
side down because the inner and the outer dipole changed 
simultaneously their direction. 

Up to this point we discussed the Poynting flux near 
the antibonding dipole resonance of the dissipationless 
model particle. Since the dielectric functions for the core 


and the shell are real (see, Eq. ([T^) only radiative damp¬ 
ing is thus accounted for so far. In a real particle, how¬ 
ever, dissipative losses occur as well. The dielectric func¬ 
tions for the core and the shell of the CaO/AUOs particle 
plotted in the left bottom panel of Fig. have imaginary 
parts. From the work of Wang and coworker^^ we know 
that the topology of the Poynting vector is very sensitive 
to dissipation. We expect therefore the Poynting flux 
of the CaO/AUOs particle to deviate from the Poynting 
flux of the model particle. Indeed for the antibonding 
resonance we found by artificially adding an imaginary 
part to the dielectric function of the model particle that 
large imaginary parts of the core and the shell dielectric 
functions lead in the shell region to a dominance of the 
inner dipole and to a strong damping of the outer dipole. 
Since the saddle points are an outgrowth of the relative 
strength of the two dipoles it is clear that not all the 
singular points of the dissipationless model particle may 
appear for a particle with dissipation. For the particular 
case of a CaO/AUOa particle this can be seen in Fig. 
where we plot for wave numbers near the antibonding 
dipole resonance in panel (a) the positions of the outer 
singular points and in panel (b) the positions of the inner 
singular points. Compared to the model particle many 
outer singular points are missing, while most of the inner 
singular points are still present. 
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Figure 9: (Color online) Top left: Position of the inside singular points (Arabic letters) in the Poynting field of the model 
particle (/ = 0.7 and a = 0.1/im) in the vicinity of the bonding resonance at « 1095.669 cm“^. Roman numbers indicate 
spectral regions with a fixed topology while the shaded area denotes the shell. The upper right and the lower left panel show 
the Poynting field for two representative wave numbers whereas the panel on the lower right is just a zoom-in into the panel 
to its left. The modulus of the Poynting vector normalized to the flux of incident energy is encoded in the color of the arrows. 
The intersections of the isoclines Seir^O) = 0 (green) and Sr{r,0) — 0 (blue) determine the singular points. The shell and the 
core are indicated by brown circles and the optical whirlpools are shown as light blue circles. 


2. Bonding dipole resonance 

Let us return to the model particle and discuss the 
topology of the Poynting vector field near the bonding 
dipole resonance at « 1095.669cm“^. The results 
are summarized in Fig. As it was the case for the 
antibonding resonance, optical whirlpools appear at the 
core-shell and the shell-vacuum interface. They corre¬ 
spond to surface localized electromagnetic modes. How¬ 


ever, since at the bonding resonance the inner and the 
outer dipole are oriented parallel, the whirlpools at the 
two interfaces rotate now in the same direction. For 
A“^ < AjC^, for instance, in region Vila, the whirlpools at 
x/a = —\^z/a = 0 and x/a = —0.1,z/a = 0 both rotate 
clockwise while the whirlpools at x/a = 0.7, z/a = 0 and 
x/a = \,z/a = 0 both rotate anti-clockwise. Behind the 
resonance, that is, for A“^ > for instance, in region 
VIIc the situation is reversed. 
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A,'^ (cm'^) (cm'^) 

Figure 10: (Color online) Position of the singular points (Arabic letters) of the Poynting field outside (a) and inside (b) 
of a CaO/AbOs particle with filling factor / = 0.7 and radius a = 0.4/rm in the vicinity of the bonding resonance at 
AJCj « 709.63 cni“^. Roman numbers indicate regions with a fixed topology which were also found for the model particle and 
the shaded regions denote the shell. 


In the vicinity of two rectified optical whirlpools energy 
flows in the opposite direction. As a result the singular 
points N and O appear in region Vila and region Vllb. 
They are caused by a S'g = 0 isocline inside the shell, 
indicated by a green line, which forms a circle around the 
core. The two singular points N and O, which separate 
the field lines entering the core from the field lines leaving 
the particle, operate thus in the same way as the two 
singular points A and B of region I outside the particle 
(see lower right panel in Fig. which separate the field 
lines entering the particle from the ones leaving it. 

In addition to the points N and O for wave numbers 
close to the resonance the saddle point D appears inside 
the particle. The closer the wave number is to the res¬ 
onance the more moves D inside the particle. In region 
Vila D is inside the shell while in region Vllb D is in¬ 
side the core. Underneath the Sr = 0 isocline at z < 0 is 
moreover hidden another Sg = 0 isocline which is coming 
from the outside. Hence, if the saddle point D is deep 
enough inside the particle an avoided crossing occurs be¬ 
tween the two Sg = 0 isoclines. In the lower left panel 
of Fig. 1^ only the Sg = 0 isocline encircling the core is 
visible. To clarify the situation we replotted it therefore 
on a magnified scale in the lower right panel. There it 
can be seen that the Sg = 0 isocline coming from the 
outside of the particle is going downwards in the shell, 
while the other Sg = 0 isocline is going towards the core. 
As a result the singular points Q and P appear. While 
the point Q is an unstable foci, that is, an energy sink 
in the ccz-plane, the point P is a saddle point. In region 
VI, where the saddle point D is still inside 

the shell but lies now above the Sg = 0 isocline existing 


in the shell and encircling the core. The topology in this 
region is the same as in region Vila, but upside down 
and with different orientation of the optical whirlpools, 
since the inner and outer dipole simultaneously changed 
directions. 

Having discussed for the dissipationless model parti¬ 
cle the topology of the inside Poynting vector field near 
the bonding dipole resonance, we now comment on the 
energy flux near such a resonance in a real particle with 
finite dissipation, taking again the Ca 0 /Al 203 particle 
as an example. For the bonding resonance we empiri¬ 
cally found by varying the imaginary parts of the dielec¬ 
tric function the outer dipole to be rather stable, even for 
large values of e” 2 - The inner dipole however disappears. 
As a result there are also no surface localized electromag¬ 
netic modes at the core-shell interface anymore. The re¬ 
sults for the bonding dipole resonance of the Ca 0 /Al 203 
particle shown in Fig. are consistent with these em¬ 
pirical findings. Whereas almost all singular points in 
the Poynting flux outside the particle are preserved be¬ 
cause of the robustness of the outer dipole, inside the 
shell most of the singularities disappeared. The situa¬ 
tion is thus just reversed to the one we found for the 
antibonding resonance. 

C. Electric field 

We now investigate for the model particle the electric 
field for wavelengths in the vicinity of the bonding and 
antibonding dipole resonance and contrast it with the 
field in the vicinity of the core resonance. For comparison 
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we also include the field near the anomalous resonance of 
an homogeneous particle. For all four resonances we cal¬ 
culate the electric field amplitude \E\'^ and the real part 
of the electric field Re(i5) for wave numbers where the 
topology of the Poynting flux is of type I shown in Fig. 
For the homogeneous particle = 1406.5cm“^, while 
for the core-shell particle = 306.026cm“^ (core res¬ 
onance), A“^ = 1095.74cm“^ (bonding resonance), and 
A“^ = 1767.375cm“^ (antibonding resonance). 


The results for the electric field are summarized in 
Fig. 11 where the value of the electric field amplitude is 
encoded in colors and the direction of the real part of the 
electric field is given by black arrows. The electric field 
of the incident wave is polarized in x-direction. Hence, 
the induced charge displacements are also in x-direction. 

Let us first discuss the field for the anomalous reso¬ 
nance of the homogeneous particle shown in panel (a) 
of Fig. [T^ From the electric field structure outside 


the particle it is clear that a dipole mode is excited. 
The field lines by definition go from positive to negative 
charges but the dipole of course points towards positive 
charges. Hence, the excited dipole is oriented in negative 
x-direction implying polarization-induced surface charges 
to be positive for the x < 0 hemisphere and negative for 
the X > 0 hemisphere. As can be seen from the linear 
scale of the plot, the variation of the electric field ampli¬ 
tude inside the homogeneous particle is very small. Yet, 
a shallow minimum occurs at the center of the particle. 

The electric field for the core-shell particle is plotted 
in panels (b) to (d) for the core, the bonding, and the 
antibonding resonance, respectively. Inside the core the 
electric field resembles for all three resonances the field of 
an homogeneous particle. On the logarithmic scale used 
for the plots this cannot be seen but a linear scale would 
reveal an onion-like structure similar to the one depicted 
in panel (a). Only for the core resonance (b) however 
would the minimum of the field still be at the center 
of the particle whereas for the bonding and antibonding 
resonance the field at the center is in fact maximal. The 
core resonance (b) has most in common with the anoma¬ 
lous resonance of the homogeneous particle (a) because 
it arises also from the resonance condition of the core di¬ 
electric function. Only a single dipole is thus excited as 
can be seen from the field lines. They originate from the 
X < 0 hemisphere of the core-shell interface and termi¬ 
nate in its X > 0 hemisphere, signalling positive and neg¬ 
ative polarization-induced surface charges, respectively. 
At the shell-vacuum interface the field lines are only re¬ 
fracted leading to an inhomogeneous field inside the shell 
which is in most parts much smaller than in the core. 

For the bonding and antibonding resonance the situ¬ 
ation is different because now two dipoles are involved. 
From the polarization-induced surface charges shown in 
Fig. [3] we already know that for the bonding resonance 
the dipoles are orientated parallel whereas for the anti¬ 
bonding resonance they are oriented antiparallel. The 
field lines shown in the panels (c) and (d) verify this. 
Let us first look at panel (c) displaying the field for the 


bonding dipole resonance. The shell-core interface in the 
X < 0 hemisphere is a source for field lines, indicating 
positive charges, whereas the shell-core interface in the 
X > 0 hemisphere is a sink for field lines, indicating nega¬ 
tive charges. At the shell-vacuum interface the situation 
is the same. Both the inner and the outer dipole are thus 
oriented in negative x-direction. Inside the shell the two 
dipole fields interfere which destroys the homogeneity of 
the field. For z = 0 and x/a « ±0.9 destructive inter¬ 
ference occurs as the outgoing field lines from the inner 
dipole cancel the inner field lines of the outer dipole. At 
zja = ±0.9 and x = 0, on the other hand, constructive 
interference occurs, because outer field lines of the inner 
dipole and inner field lines of the outer dipole run par¬ 
allel leading to a high electric field. At the antibonding 
resonance shown in panel (d) the shell-core interface in 
the X < 0 hemisphere is a sink for field lines while in 
the X > 0 hemisphere it is a source. The shell-vacuum 
interface on the other hand is a source for field lines in 
the X < 0 hemisphere and a sink x > 0 hemisphere. The 
inner and the outer dipole are thus oriented antiparallel. 
In the shell the inner field lines of the onter dipole inter¬ 
fere with the outer field lines of the inner dipole, leading 
at z/a = ±0.9 and x/a = 0 to a minimum in the field. 

The electric field in core-shell particles differs therefore 
dramatically from the one in an homogeneous particle. 
Whereas the latter is homogeneous the former is very in¬ 
homogeneous, varying in the shell for the bonding and 
antibonding resonance over almost two orders of magni¬ 
tude. At the bonding resonance the field amplitude is 
also high in the core but for the antibonding resonance 
it is one order of magnitude smaller in the core than in 
most parts of the shell. The inhomogeneity of the elec¬ 
tric field can be used to control the dissipation of energy 
inside a core-shell particle by varying the wavelength of 
the incident light. For the model particle discussed in 
this subsection there is of course no dissipation since the 
dielectric functions of the core and the shell are purely 
real but for physical particles the dielectric functions are 
complex and a control mechanism is conceivable. 


D. Dissipation 


In this section we present results for the dissipation 
of energy inside a Ca0/Al2 03 particle. According to 
Eq. (15) it depends on the imaginary parts of the dielec¬ 


tric functions as well as on the amplitude of the electric 
field. Since the latter is strongly enhanced in the vicin¬ 
ity of anomalous resonances, even small imaginary parts 
of the dielectric functions may thus lead to a huge dissi- 
pation.12^ The field enhancement in a core-shell particle 
is moreover rather inhomogeneous, there will be hence 
spots inside it where dissipation is strong and spots where 
it will be rather weak, irrespective of the magnitude of 
the (spatially homogeneous) imaginary parts of the di¬ 
electric functions. 


In Fig. 12 we show the electric field and the dissipa- 
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Figure 11: (Color online) Electric field amplitude \E\'^ and the direction of the real part of the electric field in the vicinity 
of the core (b), the bonding (c), and the antibonding (d) dipole resonance of the model particle (/ = 0.7 and a = 0.1pm). 
Panel (a) shows for comparison results for the dipole resonance of a homogeneous particle of the same size made out of the 
core material. Black arrows denote the direction of the real part of the electric field while the colors give the field amplitude. 
The wave numbers used for the plots are given in the main text. 


tion for the Ca 0 /Al 203 particle near the bonding and 
the antibonding dipole resonance. The wave numbers 
are in both cases in the topological region I where the 
Poynting flux outside the particle displays Tribelsky ears. 
For the bonding resonance = 715 cm“^ leading to 
e'l = 0.104 and 62 = 0.304, while for the antibonding 
resonance = 877 cm“^ giving rise to e" = 0.0492 
and £' 2 ' = 0.0657. 

The electric field inside the Ca 0 /Al 203 particle shown 
in the upper two panels of Fig. is inhomogeneous as 
it was for the model particle discussed in the previous 
section. Due to constructive and destructive interference 
of the electric helds of the excited dipoles inside the shell 
there are points where the amplitude of the electric field 


is enhanced and points where it is very small. For the 
bonding resonance (panel (a)) the field amplitude inside 
the core is high and more or less homogeneous whereas 
inside the shell the held amplitude varies over two orders 
of magnitude reaching at x/a = ±0.9 and z/a = 0 its 
smallest value which is two orders of magnitude lower 
than inside the core. Compared to the dissipationless 
model particle the electric held inside the core of the 
Ca 0 /Al 203 particle is no longer symmetric along the x- 
axis and it reaches its maximum no longer in the center 
of the particle but around z/a — 0.7 and x/a = 0. On 
the logarithmic scale this cannot be seen but a linear 
scale would reveal these features. The spatial distribu¬ 
tion of the dissipation near the bonding resonance (panel 
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Figure 12: (Color online) The top two panels show the electric field amplitude \E\^ and the direction of the real part of the 
electric field in the vicinity of the bonding (a) and the antibonding (b) dipole resonance of a Ca 0 /Al 203 particle with / = 0.7 
and a = The wave numbers are respectively = 715 cm“^ and = 877cm“^. The two bottom panels depict the 

dissipation inside the particle for the bonding (c) and antibonding (d) dipole resonance, respectively. Black arrows in the upper 
two panels denote the direction of the real part of the electric field while colors indicate the electric field amplitude (upper two 
panels) or the dissipation (lower two panels). 


(c)) parallels more or less the field distribution. It is ho¬ 
mogeneous inside the core, where it is also rather high, 
and varies strongly inside the shell reaching its smallest 
value at x/a = ±0.9 and z/a = 0. The results for the 
antibonding resonance, given in panels (b) and (d) of 
Fig. [T^ show the same overall trend. Inside the core the 
field and the dissipation are homogeneous whereas in the 
shell strong variations occur. In contrast to the bonding 
resonance the electric field inside the core is now however 
smaller than in most parts of the shell where destructive 
and constructive interference lead moreover to a strongly 
varying field. This feature is also preserved in the spatial 
distribution of the dissipation. It varies strongly within 


the shell and is in most parts of it two orders of magni¬ 
tude higher than inside the core. There are however also 
narrow regions where dissipation is strongly suppressed. 

The spatial inhomogeneities of the electric field in the 
shell of the model particle can thus be also found in 
the shell of the Ca 0 /Al 203 particle. Since the dielec¬ 
tric functions of the latter have finite imaginary parts, 
the field inhomogeneities lead now however also to in¬ 
homogeneities in the dissipation. Within the shell the 
dissipation of energy varies for both the bonding and the 
antibonding resonance over almost two orders of magni¬ 
tude giving rise to extended regions of high dissipation 
and well localized cold spots where dissipation is rather 
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low although the dielectric function of the shell is spa¬ 
tially homogeneous. 


IV. CONCLUSIONS 

In our previous work we proposed to use the antibond¬ 
ing dipole resonance of dielectric core-shell particles hav¬ 
ing a core with negative and a shell with positive elec¬ 
tron affinity to determine by infrared attenuation spec¬ 
troscopy the number of surplus electrons of the particle. 
Due to the particular choice of the electron affinities the 
electrons would be confined to the shell region leading 
to a high volume electron density which in turn would 
strongly affect the electric polarizability of the particle 
and hence the position of the extinction resonances. In 
particular the maximum of the antibonding dipole reso¬ 
nance, which we initially called shell resonance since it is 
absent for homogeneous particles, turned out to be very 
charge sensitive. To understand the mechanism leading 
to this resonance and also to make contact to the field of 
plasmonics, we investigated in the present work neutral 
dielectric core-shell particles, using a Ca 0 /Al 203 par¬ 
ticle as a physical example and a dissipationless model 
core-shell particle as an idealized system. The scattering 
of light is in both cases mediated by infrared transverse 
optical phonons. 

We showed that the series of bonding and antibond¬ 
ing resonances arising in the surface mode regime of the 
shell are attached to in-phase and out-of-phase multipo¬ 
lar polarization-induced surface charges building up at 
the two interfaces to satisfy the boundary conditions of 
the Maxwell equations. For the series of bonding reso¬ 
nances the surface charges are in-phase whereas for the 
series of antibonding resonances they are out-of-phase. 
The inner and outer dipoles associated with the two low¬ 
est resonances are thus aligned parallel for the bond¬ 
ing and antiparallel for the antibonding resonance. We 
demonstrated that both series of resonances scatter light 
anomalously, that is, originate from selective excitation 
of surface modes giving then rise to an inverse hierar¬ 
chy in the extinction cross section. We also analyzed the 
topology of the Poynting flux in the vicinity of the bond¬ 
ing and the antibonding dipole resonance. The outside 
Poynting flux near these two resonances is determined by 
an effective dipole describing the effective outer field of 
the inner and outer dipole. The topology of the outside 
Poynting flux is thus for the bonding and antibonding res¬ 
onance similar to the Poynting flux of an homogeneous 
particle. Inside the particle however the Poynting flux 
is determined by the interplay of the outer field of the 
inner and the inner field of the outer dipole. As a result, 
the topology of the inside Poynting flux depends on the 
alignment of the two dipoles and is thus different for the 
bonding and antibonding resonance. 

Inside the particle the Poynting flux, the electric field, 
and the dissipation of energy are highly inhomogeneous. 
Of particular interest, perhaps also from the technolog¬ 


ical view point of laser heating, is the inhomogeneous 
dissipation near the antibonding dipole resonance. In¬ 
side the shell there are broad banana-shaped hot spots 
where the dissipation is almost two orders of magnitude 
larger than in the core as well as rather localized cold 
spots where the dissipation is almost an order smaller. 
Although we do not expect the enhanced dissipation in 
the shell to be a problem for the optical charge measure¬ 
ment we proposed since the surplus electrons are bound 
by an energy on the order of the shell material’s electron 
affinity, which is for technologically relevant dielectrics 
such as AI 2 O 3 in the range of electron volts, the heating 
of the confined electron gas should be investigated in the 
future. 
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Appendix 

Usin g the orthogonality of the vector spherical har- 
monics,l2Sl2ll the eight expansion coefhcients a„, 6„, c„, 
dm fn, dn, Vn and Wn Can be straightforwardly calcu¬ 


lated from the core-shell boundary conditions at r = b, 

(A 2 -Ai)xr = 0 , ( 17 ) 

(i 72 -i?i) xr = 0 , ( 18 ) 

and the shell-vacuum boundary conditions at r = a 

(As-I-Ain - A 2 ) X r = 0 , ( 19 ) 

(i?,+i?i„-i? 2 ) xr = 0 , ( 20 ) 

yielding eight equations for the eight unknowns. 


The coefficients a„ and determining the scatter¬ 
ing components of the electromagnetic fields in the outer 
space, Eqs. turn then out to b^^ZEll 

a-n = C~^^i’n{y)Wn{'m2y) “ A„x'„(w2y)] 

- m2'4)'n{y)[i’n{m2y) - A„Xn(w2y)]| , (21) 

K = D;;;i|m2V’n(y)[V’U”^22/) - S„x'(m2y)] 

- '>P'n{y)[i’n{rn2y) - BnXn{jn2y)]\ ( 22 ) 
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with 


with 


^ _ m2'ipn{m2x)'tp'^{mxx) - nii'ip'^{m2x)iin{mix) 

m2Xn{rn2x)'ip'^{mix) - mix'n{'rn2x)'tpn{mix) 

(23) 

^ _ m2ipn{mix)'tp'^{m2x) - ■miipri{rn2x)'ip'^{mix) 

m2Xni'i^2x)i’n{rriix) - rni'ip'„{rnix)xni'm2x) 

(24) 

Cn = Uiy)['(Pl,{m2y) - A„x'„(’7^2y)] - w 2 C(y) 

X [V'n(w22/) - AnXnim2y)] , (25) 

Dn = m2^n{yMn{m2y) “ B„x'„(™ 22 /)] “ Civ) 

X [i’nim2y) - BnXnim2y)] ■ (26) 


Fn = H’niy)^n{y) “ V'n(y)?n (Z/)] . (29) 


while for the fields inside the shell, Eqs. 0-§ , the co¬ 
efficients become 


For the fields inside the core, Eqs. ([I0|)- 0 , the coeffi¬ 
cients read 


c„ = -£>„ ^Fnm 2 mi 

^ [(Pn{m2x)x'„{m2x) - Xu{m2x)ip'„{m2x)] 

[miXnim2x)ip’„{mix) - m2'4’nirnix)x’„im2x)] 

(27) 


dn = -C^ ^Fnm2mi 


[^n{m2x)Xn{m2x) - Xn{m2x)'tlj'„{m2x)] 
[m2Xn{m2x)'ip'„{mix) - mi'ipn{mix)Xn(.'n^2x)] 

(28) 


/„ = D~^m2Fn , (30) 

g„ = C~^m2Fn , (31) 

= D~^m2FnBn , (32) 

Wn = C~^m2FnAn ■ (33) 
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